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We extend our studies of thermal diffusion of non-topological solitons to anharmonic FPU-type 
chains with additional long-range couplings. The observed superdiffusive behavior in the case of 
nearest neighbor interaction (NNI) turns out to be the dominating mechanism for the soliton dif- 
fusion on chains with long-range interactions (LRI). Using a collective variable technique in the 
framework of a variational analysis for the continuum approximation of the chain, we derive a set 
of stochastic integro-differential equations for the collective variables (CV) soliton position and the 
inverse soliton width. This set can be reduced to a statistically equivalent set of Langevin-type 
equations for the CV, which shares the same Fokker-Planck equation. The solution of the Langevin 
set and the Langevin dynamics simulations of the discrete system agree well and demonstrate that 
the variance of the soliton increases stronger than linearly with time (superdiffusion) . This result 
for the soliton diffusion on anharmonic chains with long-range interactions reinforces the conjecture 
that superdiffusion is a generic feature of non-topological solitons. 



I. INTRODUCTION 

Solitons appear in various different fields of physics 
such as solid state physics, quantum theory, nonlinear op- 
tics, fluid dynamics, biophysics or even cosmology 0,0- 
The numerical experiment done by Fermi, Pasta and 
Ulam in 1955, investigating the thermalization of a one- 
dimensional chain with anharmonic interaction potential, 
was certainly one of the key events in nonlinear science 
Hereafter, FPU chains, despite of their simplicity, 
served as frequently used models for analyzing the prop- 
erties of microscopic structures, especially in computer 
simulations. More than fifty years after the above men- 
tioned experiment by Fermi, Pasta and Ulam, complete 
understanding of this putatively simple system, which 
still reveals new astonishing phenomena is yet to be 
reached. The thermalization of the FPU system, the 
heat conduction in FPU systems or the properties of lat- 
tice solitons and discrete breathers are still tasking many 
physicists (see for example CHAOS Volume 15: focus 
issue: THE " FERMI-PASTA-ULAM" PROBLEM-THE 
FIRST 50 YEARS). 

Physicists usually use simplified models to describe and 
explain the most important features of more complex sys- 
tems. An example is the "invasion" of physicists into 
biology, where they aimed to explain some of the most 
interesting processes in living cells such as the DNA dy- 
namics (DNA melting, DNA transcription) or the func- 
tionality of proteins as molecular motors, and the stor- 
age and transport of energy released by the hydroly- 
sis of adenosine-triphosphate (ATP) d B 0, S M E3 • 
Oversimplified models, which were formerly used to pre- 
dict freely moving solitons along DNA or protein a- 
helices may be too rough to explain the dynamics of 
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biomolecules. But, for example, taking the role of dis- 
creteness, thermal perturbations and the underlying ge- 
ometry of the system into account, it is possiblee to ob- 
tain good agreements with the results of experiments, 
which today work at the scale of single molecules. 
The Davydov model [ToL El for energy transport and 
storage in proteins has attracted much attention in the 
past. Davydov proposed a self-trapping of vibrational 
energy of the C=0 vibration, similar to the concept of 
polarons in solid-state physics. Recently, a direct obser- 
vation of self-trapped states in pump-probe experiments 
of a model protein (ACN) , showed that the self-trapping 
could describe some important features of proteins (al- 
though the lifetime of the excitation in isolated molecules 
is shorter than expected by Davydov) [l2|. In a differ- 
ent approach by Yomosa a molecular mechanism of 
muscle contraction was introduced, which assumes that 
the H-bonding interaction of the peptide groups (which 
stabilize the helix structure) can be modeled by a FPU- 
type lattice. The energy transport along the helix was 
described in the continuum approximation by use of the 
stability of the Korteweg-de Vries (KdV) or the Toda soli- 
ton |l4l IT5I flq | . Although the continuum approximation 
and the disregard of the helical structure seems to be bold 
to model real processes, the concept of lattice solitons in 
protein backbones was further refined. It was numerically 
demonstrated that protein backbones with a more realis- 
tic structure, support stable non-topological solitons [l7| . 
Recently, the concept of solitons along hydrogen bonds 
in a helices was applied to realistic models for proteins 
in order to explain Raman measurements, which identify 
excitations in proteins, which can not stem from vibra- 
tional modes in the molecule and are too long-lived to 
result from large-scale vibrations of the protein [Isl IT^ | . 
The three-dimensional model of the a helix consists of 
three chains with hydrogen bonds (secondary structure), 
which arc connected with peptide bonds (primary struc- 
ture), modeled by Lennard- Jones or Toda potentials. It 



2 



was demonstrated that a perturbation in one of the three 
H-bond lines can create a triple-soliton in the three hy- 
drogen bonds, which travels (phase-locked by the peptide 
bonds) along the chain and would present a long-lived en- 
ergy source for the spectroscopically observed phonons at 
energies of about ~ 100 crrT 1 . 

Often, the dynamics of more realistic models can only 
be examined or dealt with in computer simulations. Yet, 
sometimes simplified models, such as the FPU-chain, are 
able to provide proper understanding of these models. 
Therefore, it is important to clarify the effect of exten- 
sions on systems such as FPU-chain and if possible, ap- 
ply analytical calculations. In the case of biomolccules, 
especially damping and thermal fluctuations appear as 
natural perturbations. The influence of the thermal 
fluctuations on solitons on FPU chains (with cubic or 
quart ic anharmonicity) has recently been investigated 
pol I21I I22I I23I ]. It was demonstrated that both types 
of non-topological solitons, pulse solitons and envelope 
solitons, show a position variance growing stronger than 
linearly in time (superdiffusion) , whereas the mechanisms 
and the time-dependencies differ. Envelope solitons are 
generally more sensitive to damping than pulse solitons. 
This sensitivity leads to an exponentially growing po- 
sition variance, which ends when the envelope soliton 
is destroyed by damping and noise. The pulse solitons 
have a much longer lifetime and the time-dependence of 
the position variance proved to have a linear (normal) 
and a quadratic (superdiffusive) part ■ The superdif- 
fusion appears to be typical for non-topological solitons 
because, whereas it was also observed for pulse solitons 
in classical Heisenberg systems 0> UK- it is absent in 
the case of topological solitons in the sine-Gordon sys- 
tem [2^. 

In this article, the aim is to extend the results to chains 
with long-range interactions (LRI). LRI are necessary 
extensions which help make the models more realistic. 
In biomolecules, charged groups or dipole moments (for 
example the dipole moment of the base pairs in DNA) 
should have an effect on solitonic or discrete excitations 
Hi HzL HH • I n the case °f the FPU system, it was demon- 
strated that chains with Kac-Baker or power-law long- 
range interactions support pulse solitons as well |28Ll3l| . 
In the case of Kac-Baker LRI, two different velocity 
branches for soliton solutions, which are separated by a 
gap, appear. Low-velocity solutions with a soliton width 
larger than the interaction radius are pulse-shaped, but 
develop a cusp when the soliton velocity is close to the 
critical velocity at the gap 29, 30] . High- velocity solitons 
with a width in the range of the interaction radius have a 
multi-component structure, where the center of the soli- 
ton is mainly determined b y th e NNI, whereas the tails 
are determined by the LRI |3l], |32| • 
According to our research, this is the first known attempt 
to study the effect of thermal fluctuations on solitons, 
of wich the shapes are governed by long-range interac- 
tions. Our goal is to investigate the consequences of the 
long-range interactions on the stochastic behavior of the 



soliton, and specifically on the time-dependence of its 
position variance. 



II. THE MODEL 

Our model is a one-dimensional chain of equally spaced 
particles of mass m (m = 1) with interatomic spacing a 
(a = 1). The displacement of the particle n from its equi- 
librium position is denoted The potential consists 
of a part Unn, with only NNI similar to the original 
potential chosen by Fermi, Pasta and Ulam, and a long- 
range part Ulr with harmonic coupling J nm between the 
particles n and m: 



L — T — Unn — Ulr. 



(1) 



with 



T =\^[-jr) , Unn = J2v(u n+1 -u n ) 



U LR = ^^2jn,m(u n - U m ) 2 . (2) 



The potential V(r) between nearest neighbors represents 
the expansion of the Toda potential to third order 



V(r) = -r 2 - -r 3 
W 2 3 



(3) 



The long-range coupling between sites n and m is of the 
Kac-Baker form |33j,|3J], a" 1 is the interaction radius of 
the LRI 



1 N 



(4) 



For modeling the environment, we choose hydrodynam- 
ical damping instead of Stokes damping (because it was 
demonstrated in [2^ that the long-wave components of 
the Fourier spectrum are strongly damped for Stokes 
damping which leads to deformations of the soliton) and 
a Gaussian white noise term which fulfills the fluctuation- 
dissipation-theorem (D hy = 2vh y T with fog set to unity). 
The equations of motion in relative displacement coordi- 
nates w n = itn+1 — u n then take the form 

Wn = V'{w n+1 ) - 2V'(w n ) + V\w n -i) 



J 



- -(e a -l)^(w n -w n+ „ l )e 01 

m 

+ Vhy(w n +1 - 2w n + W„_i) 



|m| 



(•5) 



In order to apply similar analytic tools as used in 
[20T l2l| . we follow Ref.[3(J and derive in the quasi- 
continuum approximation (QCA) a partial differential 
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equation (PDE), which is of the Boussinesq (Bq) -type 



dtw = (c 2 1) 



k 2 — d 2 . 



w + d 2 x {w - w 2 ) + Xd 2 x diw 



-v h yd 2 x dtw + \fD^d 2 x t{x, t) (6) 



FIG. 1: (Color online) Soliton solutions of eq. © with the 
same normalized velocity c D = v/c = 1.01 for a chain with 
nearest neighbor interactions (J = 0) and for a chain with 
a relatively weak long-range coupling J = 0.1 for different 
values of the interaction radius a -1 . 
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(7) 
(8) 



with an additional damping, noise and long-range term 
(a pseudo differential operator) which is defined as 



(« 2 



1 f°° 

dl) ~ lf{x) =2^j d * ,e ~ alX ~ X] f( x ') ■ ( 9 ) 



In [3fJ, |31| , the solution of equation JBJ (for Uh y = 
and D hy = 0) was given implicitly. It was demonstrated 
that the chain with Kac-Baker LRI supports soliton so- 
lutions in the velocity interval c < v < v c , where c is the 
velocity of sound and v c ~ v/ (4c 2 — l)/3 is the critical 
velocity, where the solitons develop a cusp like peakons, 
which are known from the theory of shallow water waves 
[37| and the integrable Camassa-Holm equation |2^. A 
stability analysis in [3(j proved that the soliton becomes 
unstable before it adopts peakon shape. It was further 
demonstrated that it is possible to construct long-range 
potentials in nonlinear lattices which support (discrete) 
peakon solutions |3^ |. 

Here, the aim is not to directly solve equation ©, but 
rather to benefit from the fact that as long as the soliton 
velocity is smaller than v c , the soliton can be relatively 
well approximated by a pulse-shaped collective variable 
(CV) ansatz. As demonstrated in Fig. ^ an increasing 
long-range radius a -1 yields larger and broader excita- 
tions even when the long-range coupling J is relatively 
small. For a = 0.5 and J = 0.1, the soliton with the nor- 
malized velocity c c = v/c = 1.01 is much broader than 
the soliton solution without long-range forces (J = 0) 
with the same c . Therefore, the continuum approxima- 
tion seems to be a very appropriate analytical tool to 
investigate long-range effects. 

The case of nearest neighbor interactions is covered by 
our theory in the limit J — > (for a^0). Here, the long- 
range term in equation © disappears when (c 2 — 1) — > 
and we obtain an improved Boussinesq equation (IBq) 
with c = 1 |2^|. The limit a — > oo leads to a IBq with 
an additional nearest neighbor coupling J/2 and a sound 
velocity c = \/l + J/2. 



III. COLLECTIVE VARIABLES 

A collective variable theory depends crucially on the 
trial function one chooses to derive ordinary differen- 
tial equations (ODEs) for the soliton parameters, such 



O—B J=0,C=1 

1=0.1, a=2. c=1.037 
— o J=0.1, a=l, c=1.082 
1=0.1, a=0.75, c=1.125 
J=0.1,a=0.5, c=1.232 




FIG. 2: (Color online) Soliton solution for different values of 
c for J = 0.1 and a = 0.3. For c = v c /c — 1.092, the soliton 
solution takes the form of a peakon. 




as width and position. As seen in Fig. [5J the shape of 
the soliton is rather pulse-like up to values c a < 1.07. A 
CV ansatz of this form has already achieved qualitatively 
good results with simulations for a large range of possible 
soliton velocities [2^. In absolute displacement coordi- 
nates u(x, t), the pulse-shaped soliton solution (in rela- 
tive displacement coordinates w(x, t) = d x u{x,t)) takes 
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the form of a kink, where the amplitude 

2A Q = / w(x,t)dx (10) 



is a conserved quantity (such as for the Bq equation). 
Therefore, we choose a trial function of the form 

u(x, t) = A tanh( 1 (t)(x - X(t))) (11) 
w(x,t) = d x u{x,t) = A ~/{t)sech 2 (~f(t)(x - X{t)) (12) 

with the time-dependent collective coordinates 7(f) (in- 
verse width) and X(t) (position) and the kink amplitude 
2A of the unperturbed initial soliton. This CV ansatz 
yields the correct results (the Bq soliton solution) in the 
case of an unperturbed system with only NNI 



V3 



r \/v 1 — 1 v •) — 

7 = v3 , A a = — z - lv 

v 2 

X(t) = vt + X a 



(13) 



In order to derive ODEs for the collective variables "f(i) 
and X(€), we apply a variational technique similar to the 
method of Whitham [371 • It is well known that a La- 
grange density of Boussinesq-type equations only exists 
for the absolute coordinate field u(x, t) 41]. In order 
to include the damping term, the generalized Hamilton 
principle of Ostrovsky et al. |40l |41| , who extended the 
method of Whitham to nonconservative systems, is here 
applied. The equation © in absolute displacement co- 
ordinates 



u - d 2 x u + 2(d x u)(d 2 x u) 
VD^d x £(x,t) = v hy dld t u . (14) 



can be written as 

dC d dC d dC d 2 dC 



du dt du t dx du x dxd t du xt 



$ (15) 



with 



it? u 2 ul A 



r = y- Y + y + \<t + Vd*u*Z(x, t) 



(c 2 - 1)k 



e- Klx - s \u(x,t)d 2 u(s,t)ds (16) 
$ = u hy d 2 x d t u . (17) 



Therefore, the generalized Hamiltonian principle for wave 
motion in nonlinear and dissipative media in the frame- 
work of our CV theory, is used here. Accordingly, the 
equations of motions can be derived from the averaged 
Lagrangian density and damping term by varying with 
respect to the CV (X x = ~ f (t) > X 2 = X(t)) 



5<£> d<£> d d<£> 



5Xi 



dX t 



dt d±i 



(18) 



The brackets signify a spatial integration over x. In [20j, 
where the soliton equation was the nonlinear Schrodinger 
equation, it was possible to include the damping in the 
Lagrangian density and, thus, no damping term <I> was 
needed. Further, the treatment of the damping for the 
KdV [2lJ , the treatment of the damping proved unprob- 
lematic, because a projection method was used for the 
CV [42| (this method is, however, rather extensive for 
the non-local Bq equation ©). 



IV. LANGEVIN EQUATIONS 

Substituting the ansatz (|11|) in the Lagrangian den- 
sity (|16fl and applying the variation l|18ll . results in the 
following set of stochastic integro-differential equations, 
where small terms , ~ 7 2 , ~ X and ~ 7 were neglected 
to achieve analytically manageable equations {9 = ft/27, 
x = x-X{t)) 



4 • 
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2\fD K y-i 2 



A„ 



Sech 2 [yx] Tank [72;] £(x, t)dx 
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VhyT 



(19) 



16 



2(V 2 -l) + ^A o7 -4(c 2 -l) 



6 + 6* + 9 4 ^"{6) 



2\AD^7 
A 



2^/Wv 



dxSech 2 [yx] Tank [yx] x£ (x , t) 
dxSech 2 [yx]£{x,t) = . (20) 



The values for higher time derivatives of 7(f) and X can 
be estimated from the simulation results. A classification 
of the contributions in orders of a small parameter e ~ j a 
makes obvious that these terms contribute only negligi- 
bly. It was additionally checked, that the numerical so- 
lution of system (|20|l 43] without noise yields practically 
the same results as the complete set without neglecting 
the small terms. At this point, it should be mentioned 
that the neglect of X is justified by very different argu- 
ments compared to those for systems of point particles in 
the overdamped limit, for instance in ratchet systems. In 
the following calculations, it will be demonstrated that 
the dynamics of the system depend on the reduced time 
scale t r — y 2 Vhyt- Therefore, time derivatives appear in 
an order C(e 4 ) smaller because we choose vty in the or- 
der 0(e 2 ). 

The section in the square brackets of l|20|) results from the 
long-range interactions, where VP" is the second deriva- 
tive of the digamma function. In order to get a first-order 
equation in X(t), we rewrite equation 120fl to 



1 / 2\l ' D h v f°° 
X 2 = vj ( 1 + — - — / dxSech 2 ['yx]Z(x, t) 

V d V A o J -00 



2y/D h V"f 



A 



o J — 00 



dxSech 2 [yx]Tanh[yx]x£,(x, t) I j (21) 
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with 



1 + -A ol {t) - 2(c 2 - 1) 



l Vf"(0) 



(22) 



For zero temperature, «d is the velocity of the damped 
soliton where the damping effects enter in the form of 
the soliton width j(t). To investigate the influence of 
the noise, the equations for the CV need to be written 
as Langevin-type equations. Thus, we approximate X in 
(|19|l and Vd in equation l|21|l with the initial velocity of the 
soliton v, extract the square root of l|21l) and approximate 
the r.h.s. by expanding the root to the first order in 
the small parameter V D hy . The result is a system of 
stochastic integro-differential equations for the CV: 



A 2 



with 



At 



-0.8v hy j 3 , A 2 =v d 



(24) 

Q / n/iy-,2 

Bn = -—. -Sech 2 [jx]Tanh[>yx] (25) 



B, 



2A a v 



JJhy 

B 21 =-——Sech 2 [ 7 x] (26) 

ZAqV 

Sech 2 [jx] Tanh [72;] x . (27) 



As in |2pJ, |21| , we proceed to find a statistically equiva- 
lent Langevin system (yielding the same Fokker-Planck 
equation as l|23l) in the Stratonovich interpretation [44j|L 
which is more convenient for further numerical and ana- 
lytical studies. The Langevin-system with two indepen- 
dent Gaussian white noise processes reads: 



a 2 



hi \ / e 
622 M £ 2 



(28) 



with 



ai = Ai 



1 D hy 



AAlv 2 



7 , a 2 =v d 



6ll = V5^^ 7, 



1 n 2 ^fjy^ 
b 22 = \/(r + r^r) 



6 180' ,4 



(29) 



At this point, it is possible to compare the Langevin sys- 
tem here presented, with the results only for NNI in ref- 
erence |21|. Despite the facts that the soliton diffusion 
in [23 was calculated for low-velocity solitons with a CV 
ansatz representing a KdV soliton, and that a different 
CV procedure was used, interesting similarities can be 
observed. The stochastic equation for 7(i) in (|28|l demon- 
strates, as in pl| . a damping induced broadening of the 



soliton where the drift term is proportional to the damp- 
ing constant and to the third power of 7 (7 ~ vy lv r f'\ 
plus a small correction which is proportional to D hy . 
The effect of the damping, namely the soliton broaden- 
ing mechanism, seems to be the same for solitons of the 
KdV and the Bq type. However, and more importantly, 
the long-range interactions do not seem to change this 
behavior (yet, of course they change the initial condition 
7o = l(t = 0) of the soliton). 

In Ref. jU, it is argued that the observed superdiffu- 
sive behavior is induced by the stochasticity of "f(t) and 
its effect on the soliton velocity. The fluctuations of the 
soliton shape lead to fluctuations of the soliton velocity, 
and therefore, to the superdiffusive terms in the position 
variance. For envelope solitons, however, it was shown 
that the superdiffusive mechanism differs. The reason 
is that the rapid broadening of the envelope due to the 
damping, and not the changes due to fluctuations, is the 
main contribution |20|. 

If it is assumed that, as for envelope solitons [2(|, the 
fluctuations of ^(t) are negligible and the main contribu- 
tion to the soliton diffusion results from the noise acting 
directly on the soliton center X(t), bi% = can be set 
and we can use the analytical result for the soliton width, 
which depends on the reduced time t r = j^Vhyt 



(30) 



to calculate the position variance 



Var[X] bll=0 = I b 2 22 (t')dt' 

t 

I- 1 ' 



7T 

f80 



D hy 

Alv 2 



{ {t')dt' 
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k B T5a + ^)/ 1 f 



A 2 v 2 



G 



7 (t) 3 



7o 3 



(31) 



Because higher order terms do not appear until t r < 
1, this result yields an approximately linear time- 
dependence for typical simulation times. Therefore, 
yar[V]b 11= o will be denoted as the normal diffusion of 
the soliton below. The result Var[X]b 22= o which we ob- 
tain by solving numerically the Langevin system l|28|l 
with the constraint 622 = bears only the stochastic con- 
tribution of "f(t) and can be regarded as the anomalous 
or superdiffusive part of the soliton position. Although 
the equality VarLY] = VarLY]b 11= o + Var[X]t 22= o does 
not hold, these two quantities can be used to illustrate 
the weight of normal and superdiffusive effects. 



V. SIMULATIONS 

The time integration of the equations of motion Q was 
carried out by use of the Heun method, which is widely 
used in the solution of partial differential equations or 
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difference-differential equations, coupled to either an ad- 
ditive or a multiplicative noise term |45| . Periodic bound- 
ary conditions were used to enable the running of long 
simulation times on systems of typically 2000 particles, 
clearly larger than the cut-off value for the LRI (100/a). 
For a system with LRI, the simulation times increase 
with the square of the system size N. To measure the 
diffusion of the solitons, over 100 different realizations of 
the system, which were carried out on a cluster with 48 
processors, were averaged. 

The results for the position variance Var[X(i)] mea- 
sured in the simulations and the numerical solution of 
the Langevin system (I28|) for a = 0.3 (Fig. [3J and 
a = 0.2 (Fig. 30] will be compared below. The purely 
normal and superdiffusive contributions Var[X]b 11= Q and 
Var[X]b 22 =o are also displayed to demonstrate that for 
the shown velocity range (c = 1.03, c Q = 1.05, c Q = 
1.07), the soliton diffusion Var[X] is practically equal 
to the superdiffusive part Var[X]b 22= o. The long-range 
coupling J and the inverse long-range radius a are chosen 
such as to yield a constant value for the sound velocity 
c= 1.51516. 

The results from both the simulations and Langevin 
system for a = 0.3 (Fig. |3J) and especially a — 0.2 (Fig. 
0}, generally agree well. The normal diffusion of the soli- 
ton is negligible for solitons with velocities c > 1.03. 
The smallness of the term Var[X(t)]b 11= Q shows that the 
direct action of the noise on the soliton position X(t) is 
practically negligible. 

The superdiffusive mechanism Var[X(t)]i, 22= o de- 
scribes the soliton diffusion very well and the agreement 
increasingly betters for larger values of the interaction 
radius or 1 . The fluctuations of 7(i) yield a fluctuating 
drift term Vd("f(t)) and thus a fluctuating X(t) with a 
stronger than linear variance Var[X(t)]i, 22 =o. Since the 
agreement for a = 0.2 (Fig. is even better than for 
a = 0.3 (Fig. [5J, the influence of the stochastic quan- 
tity "f(t) on the soliton velocity Vd (where the long-range 
terms contribute) can be regarded as the driving force of 
the superdiffusive mechanism. 

The Kac-Baker form keeps the sum over the coupling 
constants J m constant Therefore, the almost perfect 
agreement between Var[X]b 22= o, and the simulation re- 
sults Var[X(t)] for smaller values of a, can be regarded 
as an effect of the stronger long-range character of the 
system. The normal diffusion 1/arLY] b 11= o for the values 
of a and c a in figures |3] and 0] is only of importance for 
very small times and yields a clearly smaller slope than 
the simulation results. This discrepancy appears in pl| 
for solitons on chains with NNI as well, and it was tackled 
in [2^ by generalizing the results for the soliton diffusion 
on the Toda lattice to other anharmonic chains. It was 
demonstrated that the phonon-induced diffusion of the 
Toda soliton (when it is expressed in terms of the soliton 
characteristics (velocity, amplitude, width) of solitons on 
FPU-like chains) can be used to explain the discrepancy 
for the prediction of the normal diffusion. 
In figures [3] and 01 the dominance of the superdiffusion 



FIG. 3: Soliton diffusion for different solitons (c = 1.03, 
c = 1.05 and c = 1.07) on a chain with J = 0.1, a = 
0.3, v hy = 0.01 and T = 0.0001. The results Var[X] bll=0 
and Var[X]b 22 =o are obtained by solving the system (1281 for 
611 = and 622 = 0. The result for l / ar[X]i, 22= o illustrates 
the dominance of the superdiffusive mechanism. 
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2000 4000 6000 8000 10000 



(a) 



1=0.1, 0=0.3, c o =1.05, v=0.01,T=0.0001 

1.2 F 1 I 1 I 1 I l— 




2000 4000 6000 8000 10000 



(b) 



1=0.1, a=0.3, c o =1.07, v=0.01,T=0.0001 

1.4 F 1 I 1 I 1 I •— 




2000 4000 6000 8000 10000 



(c) 



7 



FIG. 4: Soliton diffusion for different solitons (c = 1.03, 
c = 1.05 and c = 1.07) on a chain with J = 0.0468, a — 
0.2, u hy = 0.01 and T = 0.0001. The results Var\X] bll =o 
and Var[X]b 22 -o are obtained by solving the system 1281 for 
fen = and 622 = 0. The result for Var[X] D22 =o illustrates 
the dominance of the superdiffusive mechanism. The small 
off-set in (a) results from the inaccuracy of measuring the 
soliton position for very broad solitons. 



J=0.046819, ct=0.2, v=0.01, c =1.03. T=0.0001 




(a) 

J=0.046819, c^0.2, v=0.01, c =1.05. T=0.0001 




t 



(b) 

J=0.046819, Cfc0.2, v=0.01, c =1.07. T=0.0001 




2000 4000 6000 8000 10000 



sets in earlier for larger values of c . This result can be 
qualitatively explained by the fact that the soliton dy- 
namics depends on the time scale t r l|30jl ■ which is pro- 
portional to the damping and the square of the inverse 
width. For large values of c , the width is small and the 
reduced time t r increases faster. Interestingly enough, 
the time coordinate r in the analytical result for the po- 
sition variance of the KdV solitons [2lJ appears with a 
pre-factor A which is also proportional to the damping 
and the square of the inverse width. The dependence 
of the soliton diffusion on such a time scale seems very 
typical and is ultimately a consequence of the universal 
behavior of the soliton width in the presence of damping 

To quantitatively see the influence of the long-range ra- 
dius a -1 and the time scale t r on the soliton diffusion, 
an analytical result for the position variance would have 
to be deduced. However, because the dependence of Vd 
on 9 involves a polygamma function 122(1 , this is a rather 
hard task. According to our observation, the superdif- 
fusive part of the simulation results can always be very 
well approximated by a mainly quadratic and cubic time 
dependence. This observation agrees with the result of a 
perturbation theory in |2l| that the superdiffusion yields 
a term ~ t 2 for the position variance of the KdV soliton. 
For solitons with relatively low velocities, the normal dif- 
fusion is more pronounced than for their counterparts on 
NNI chains for simulation times t < 10000. The LRI soli- 
tons become very broad for small velocities and appear 
more as a collective excitation of the system, than local- 
ized energy pulses (see FigQ]for c = 1.01). The superdif- 
fusive character for solitons, with for example c D = 1.01, 
appears only for t > 10000, when t r reaches larger values. 

The theory of section IV holds for chains with only 
nearest neighbor interactions (by setting J = 0) as well. 
In this case, the results agree well to those of earlier stud- 
ies The diffusive mechanism is basically the same 
as for LRI, but because the discreteness effects are much 
stronger, the solitons for the nearest neighbor case are 
rather narrow and can not manage velocities as high as 
for LRI. Therefore, the normal diffusion is dominant for 
broad solitons at quite small soliton velocities, for exam- 
ple c = 1.005. Superdiffusive effects can hardly be seen, 
because these low-energy solitons are rapidly destroyed 
by the noise and damping. For larger values of c D , such 
as c = 1.01 and c = 1.02 (in Figures |SJ a and b), the 
superdiffusive contribution dominates. 
The numerical results for the soliton diffusion on chains 
with Kac-Baker long-range interactions generally suggest 
that the superdiffusive mechanism is more pronounced 
than it is for nearest neighbor interactions. For LRI, 
the range of possible velocities where the superdiffusion 
dominates is clearly enlarged. The solitons in the LRI 
case also have higher energies and longer life times which 
because of the quadratic and cubic terms in the time- 
dependence of Var[X(t)], lead to the dominance of the 
superdiffusion. 



(c) 
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FIG. 5: Soliton diffusion for different soliton velocities (c = 
1.01 and Co = 1.02) on a chain with only nearest neighbor 
interactions ( J = 0, c = 1) and Vhy = 0.003 and T = 5 • 10~ 6 . 
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VI. 



(b) 



SMALL-NOISE EXPANSION AND 
LONG-TIME LIMIT 



In order to gain some deeper insight in the diffu- 
sion mechanism of the long-range solitons and the time- 
dependece of yar[X(i)], we develop a small-noise ex- 
pansion of the Langevin set (|28[1 . The most impor- 
tant step in the analytical treatment is the rather com- 
plicated dependence of the velocity Vd on the inverse 
width j(t). An expansion of the polygamma function 
^-J.-^-^ + ^for large = k/2j 
|50j is already rather appropriate for typical values of 
9{t = 0) S [1-5,3] for the considered long-range solitons 
with c = 1.51516 in Fig. [3] and Fig. This approxi- 
mation is even for smaller values of 9 justified after some 



time t > l/j 2 Vh y , because of the damping-induced grow- 
ing of the soliton width 7 -1 . It can be shown that the 
soliton velocity Vd depends in this approximation only on 
an additional term which is proportional to 7 2 : 



= Jl + l4 o7 (*)-2( C 2-l) 



02 + e^"{e) 



1 + -A ol {t) - 2(c 2 



1)[ 



6(9 2 - 



A 



hit) 



4A Q 



4(c 2 
8 Al 



1) 



2(c 2 -l), 



lit? . (32) 



The small-noise expansion of the Langevin set (|28(l fol- 
lows reference |5lJ and was already successfully applied 
to the problem of the diffusion of low-velocity solitons 
in the FPU system without long-range interactions [2l| . 
We seek an asymptotic solution of the form 



lit) 
X(t) = 



T-it) 



r co) 



-e 7 (1) (t) 



(33) 
(34) 



where e is a small parameter which is formally introduced 
to consider the influence of the noise terms as small per- 
turbations (V D hy ~ e). The contribution of the drift 
term in the different orders of e are calculated following 
the rule 



Oi(7(*)) = ad7 (0) + 



E 

m— 1 



c .m (t) daii^jt)) 
7 [) d 7 (°)(t) 



With C12 = Vd, the equations in the order e° read: 

d 7 (°)(t) = -0.8u hy ^°\t) 3 dt 

dX^(t) = (c+t^it) 

2 (c 2 - 1) 



(35) 



(36) 



81 c 3 



7(*r H« (37) 



and describe the damped soliton. The first corrections 
due to the noise terms appear in the next order e : 

d 7 (1) W = -2.4^ 7 ( 1 )(t)7 (0) (<) 2 * 
5 A a v 



-dWi 



(38) 



,/.v" ( /. = -"V)(c+i£ 
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If 



substitute the result for and 



7' 



7o/yX6i^y7p + 1 into P9"j) . we can calculate the first- 
order expression for the variance of the soliton position 
due to the noise 



VariX^t)} 



lim < Xi{t)Xx(8) > 

s— vi£ 



(40) 



This is the first order result for the soliton position vari- 
ance Var[X(t)] when we finally set e = 1 which is justi- 
fied if the noise in the system is sufficiently small. Since 
the small-noise expansion yields a very lengthy result for 
the position variance, which depends on the time scale 
t r = l.&Vhy~i1t, we prefer to demonstrate the good agree- 
ment just for one typical example (Fig. [HJl. 



FIG. 6: Simulation results for the position variance of a soli- 
ton with Co = 1.03 on a chain with J = 0.046819, a = 0.2, 
c = 1.51516, v hy = 0.01 and T = 0.0001. The results agree 
well with the numerical solution of the Langevin system 128ft . 
The result from the small-noise expansion reproduces the nu- 
merical solution of the Langevin system qualitatively well. 



1=0.46819. a=0.2, v=0.01, c =1.03, T=0.000i 



simulation 
-- Var[X(t)l 




(a) 



Further comparisons reveal that the analytical result 
for VarLY(i)] begins to fail for high velocities v close to 
the critical velocity v c . The reason is the restriction to 
the first order corrections. For high- velocity solitons, the 
width is quite narrow and the fluctuations of the width 
are strong and can no longer be well described by ^'(t). 
The long-time limit for the soliton diffusion result in the 
small-noise expansion predicts an asymtotic soliton diffu- 
sion where the underlying soliton shape determines only 
the numerical constant in 



Var[X(t)] inf 



1 D h « 

64 c 2 v 2 v]^ 



3/2 



(41) 



We can not prove this prediction in our simulations be- 
cause the result Var[X(t)] for the small-noise expansion 
reaches the asymptotic t 3 / 2 time-dependence after very 
long times. It would take times larger than 10 6 for the 



soliton of Fig. |B| For these times, the soliton is not a lo- 
calized energy pulse but smeared out over the system and 
the soliton velocity v approaches the velocity of sound c. 
Therefore, we do not assume the long-time limit to play 
a role in energy transfer mechanisms. As long as the soli- 
ton is not too broad, the time dependence of the soliton 
diffusion can be well approximated by a quadratic and 
cubic term in t. 



VII. CONCLUSIONS 

It was observed that the diffusion of solitons on chains 
with Kac-Baker long-range interactions is in many re- 
spects very similar to their counterparts on chains with 
nearest neighbour interactions. A theory which describes 
both cases very well, was derived here. Especially the 
dependence of j(t) on t r , and the superdiffusive mecha- 
nism with the mainly quadratic contribution in time to 
the position variance, in both cases, gives rise to the con- 
jecture that superdiffusion is a generic feature for non- 
topological solitons. 

In the case of long-range interactions, the solitons can 
adopt rather high velocities and it was observed that the 
superdiffusive mechanism of the Langevin system, which 
was derived in the framework of the applied CV the- 
ory for the non-local Boussinesq equation with noise and 
damping terms, describe them well. The prediction of the 
derived Langevin set for the superdiffusion of the soliton 
is especially accurate for small values of a, where soli- 
tons (for a certain velocity c ) possess more energy than 
for larger values of a. It is not surprising that the su- 
perdiffusion dominates for massive solitons because the 
direct influence of the noise on the soliton position is less 
important. The superdiffusion is more present for soli- 
tons with LRI, because the possible velocity range and 
the lifetime of the solitons are distinctly enlarged. For 
a chain with nearest neighbor coupling, solitons could 
only be described to rather moderate velocities, where 
the normal diffusion is more present. 
In the framework of a small-noise expansion for the 
Langevin system (128(1 , we derived an analytical result for 
the position variance Var[X(t)] which agreed well with 
the simulations when the soliton velocities are not too 
large. This result depends in a quite complicated way 
on the time scale t r . But the result confirms the former 
observation that the superdiffusion yields contributions 
in t 2 which turn into a t 3 / 2 dependence for long times 
when the soliton is rather broad with a velocity v close 
to the sound velocity c. 

For small times, the soliton with LRI showed a normal 
diffusion which is larger than the result Var[X (t)]b 11= o 
of the CV theory. The problem in determining the nor- 
mal diffusion of the soliton may be two-fold. In the case 
of a chain with nearest neighbor coupling 0, [2^ , it was 
established that the phonons have a significant contribu- 
tion to the diffusion of solitons, especially for low- velocity 
solitons. This result could possibly be similar for LRI 
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solitons. When, as in this paper, higher velocities are 
examined, the KdV equation is no longer suitable de- 
scribing the situation well. Instead, the Bq equation, 
where a second-order time derivative of the wave func- 
tion appears, which results in ODEs for the CVs where 
the soliton position appears quadratic {X 2 ), has to be 
used. To approximate this result by a Langevin-type 
equation, the square root of the r.h.s. of l|21|) had to be 
expanded because this procedure yields the correct veloc- 
ity v c i without noise. However, with this procedure and 
by substituting with a constant value v, the result has 



the tendency to underestimate the influence of the noise 
on the CV X. Thus, it is difficult to decide wether it 
is the phonons or the applied approximations, which are 
responsible for the observed discrepancies. 
Nevertheless, the predicted superdiffusive behavior of lo- 
calized energy pulses in systems with long-range forces 
can be observed for long enough times. Moreover, the 
long-range interactions lead to much larger soliton ener- 
gies and lifetimes, so that the soliton diffusion could be 
an observable effect in real systems. 
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